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Abstract 

In an unreliable single-hop broadcast network setting, we investigate the throughput and decoding- 
delay performance of random linear network coding as a function of the coding window size and the 
network size. Our model consists of a source transmitting packets of a single flow to a set of n users 
over independent erasure channels. The source performs random linear network coding (RLNC) over k 
(coding window size) packets and broadcasts them to the users. We note that the broadcast throughput 
of RLNC must vanish with increasing n, for any fixed k. Hence, in contrast to other works in the 
literature, we investigate how the coding window size k must scale for increasing n. Our analysis reveals 
that the coding window size of 0(ln(n)) represents a phase transition rate, below which the throughput 
converges to zero, and above which it converges to the broadcast capacity. Further, we characterize 
the asymptotic distribution of decoding delay and provide approximate expressions for the mean and 
variance of decoding delay for the scaling regime of k = w(ln(n)). These asymptotic expressions reveal 
the impact of channel correlations on the throughput and delay performance of RLNC. We also show how 
our analysis can be extended to other rateless block coding schemes such as the LT codes. Finally, we 
comment on the extension of our results to the cases of dependent channels across users and asymmetric 
channel model. 

Index Terms 

Broadcast, Delay Analysis, Erasure Channel, Network Coding. 

I. Introduction 

We consider an important transmission scenario, occurring in many communication systems, whereby 
a source must broadcast common information to many users over wireless channels in a timely manner. 
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Such a scenario occurs, for example, in a satellite or cellular network where a satellite or base station 
broadcasts a large file or streaming multi-media data to many users within their footprint over unreliable 
channels. Another example occurs in a multi-hop wireless network where each node broadcasts control 
information to all its immediate neighbors to coordinate medium access, power control, and routing 
operations. We note that such local sharing of control information (such as queue-length or other pricing 
information) is common to many provably efficient network controllers (e.g. CQ, EL etc.). 

In this work, the essential components of such wireless broadcast systems are modeled through a 
transmitter broadcasting consecutive blocks of k data packets to n users over independent and identically 
fading time-correlated erasure channels with steady state erasure probability p. Assuming that the trans- 
mitter is infinitely backlogged, we consider transmission strategies that transfer the data in blocks of k 
packets, which include the class of block coding strategies. Among all such block transmission strategies, 
it has previously been shown (see [4]) that, for any fixed n and k, the Random Linear Network Coding 
(RLNC) strategy (see Section [III] for a detailed description) asymptoticall) 3 minimizes the number of 
transmissions required to complete the transfer of all k packets at all n users (also called the block 
decoding delay). 

With this motivation, we focus on the scaling performance of RLNC as a function of k and n with 
respect to the following two key metrics: the (broadcast) throughput, defined as the data transfer rate to 
all users; and the (broadcast) decoding delay, defined as the amount of time spent between the start of 
a block transmission and its completion (i.e., successful decoding) at all the users. 

It is not difficult to see that the (broadcast) capacity of such a collection of n erasure channels, for 
any n, is equal to (1 — p) packets per time slot. Moreover, this maximum limit on the throughput can 
be arbitrarily closely achieved by encoding information into an arbitrarily large block size, k. Yet, this 
is not attractive since it leads to a decoding delay that diverges to infinity. In this work, we address 
the question of whether RLNC can achieve throughput arbitrarily close to the capacity while yielding 
acceptable decoding delay. The main contributions of this work are: 

• We find that the broadcast throughput of RLNC must vanish for any fixed k as n tends to infinity. 
We expose the cause of this behavior through a key example (see Section [TT]) , which motivates our 
search of a proper scaling of the block size k with increasing n. 

• Using upper and lower bounds on the throughput and decoding-delay performance of the broadcast 
system, we show that a phase transition in the performance of our system occurs at the block length 



'This asymptotic is with respect to increasing field size over which the data packets are defined (see Section Hill). 
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scaling ratd_| of k = ©(ln(n)) with respect to the network size n. Specifically, we show that if 
k increases slower than ln(n), then the broadcast throughput of RLNC converges to zero, and if 
k increases faster than ln(n), then the broadcast throughput of RLNC converges to the broadcast 
capacity of (1 — p). In Section [TTJ we shall note that the nature of this phase-transition is different 
from the previously observed phase-transition phenomenon in [5] due to a key difference in the 
scaling of k and the metric of focus. 

• We characterize the asymptotic distribution of decoding delay and compute the mean and variance 
of decoding delay for the scaling regime of k = w(ln(n)) using extreme value theory. 

• We provide numerical results to substantiate our findings. Our results verify the phase transition with 
respect to the scaling rate and, not surprisingly, indicate that the rate of convergence to capacity 
can be increased by choosing a faster scaling of k with respect to n at the cost of higher decoding 
delay. Our results also show that the approximate expression for mean decoding delay obtained using 
extreme value theory is accurate even for small values of n. 

These results collectively imply that RLNC can achieve throughput-delay tradeoff of ((1— p), Q(ln(n))). 
This is an attractive result as it indicates that as long as the coding block size scales super-logarithmically 
(i.e., very slowly) with the network size, the broadcast capacity is achievable with a simple policy such 
as RLNC. 

The rest of the paper is organized as follows. In Section [TTJ we overview some of the relevant work 
in this context and provide an example that motivates this work. After introducing the main system 
components in Section [TTll we provide our throughput and delay analysis of RLNC for the case of time 
invariant erasure channels in Section |W] We present some numerical results to substantiate our findings in 
the Section [V] In Section [V]J we comment on three important extensions of our results - analysis of other 
rateless block coding schemes such as the LT codes, dependent channels across users, and asymmetric 
channel model. Finally, our conclusions are provided in Section IVIIi 

II. Related Work and Motivating Example 

Our model is similar to that considered in (4J, Q, O. In 0, the authors quantify the reliability 
gain of RLNC for a fixed coding window size and show that RLNC significantly reduces the number 
of retransmissions in lossy networks compared to an end-to-end ARQ scheme. The delay performance 

2 We use the standard order notation: g(n) = o(/(n)) implies lim n ^ cxi (g(n)/ f(n)) — 0; and g(n) = u>(f(n)) implies 
lim n _ >00 (<;(n)//(n)) = oo; and g(n) = f2(/(n)) implies lim n ^oo (g(n)/f(n)) > c for some constant c; and g(n) = 0(/(n)) 
implies ci < \im n ^ 00 (g(n) / f(n)) < C2 for some constants ci and c-z- 
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gains of RLNC were observed in [4]. They show that, for a fixed coding window size k, the network 
coding capability can lead to arbitrarily better delay performance as the system parameters (number of 
users) scale when compared to traditional transmission strategies without coding. 

Also, in a similar setup as in this paper, it has recently been shown in Q that, for any given coding 
window size, there exists a phase transition with respect to decoding delay such that there is a threshold on 
the number of transmissions below which the probability that a block of coded packets can be recovered 
by all the nodes in the network is close to zero. On the other hand, if the number of transmissions is 
slightly greater than the threshold, then the probability that every node in the network is able to reconstruct 
the block quickly approaches one. 

All of the aforementioned works S, Q, study the gains of network coding as the system size 
grows while the coding window size is held constant. In particular, they show that the decoding delay 
of RLNC scales as 0(m(n)) for a fixed coding window size as n — > oo. However it can be seen that 
when the coding window size is held constant, the throughput of the system goes to zero as the system 
becomes large because each user receives a block of k packets in 0(m(n)) time slots. Therefore, it is 
important to study the system when k is scaled as a function of n. In this paper, we investigate the 
decoding delay when the coding window size scales as a function of the network size. We observe that 
there is a phase transition with respect to broadcast throughput such that when the coding window size k 
scales asymptotically slower than ln(n), the broadcast throughput converges to zero. On the other hand, 
scaling k asymptotically faster than ln(n) leads to the broadcast throughput approaching (1 — p). This 
phase transition is fundamentally different from that observed in |5] which relates to the decoding delay 
when k is fixed while the phase transition observed in our work relates to the nature of scaling of k as 
a function of network size n. 

The following example motivates our investigation of the throughput-decoding delay tradeoff of RLNC: 

Example 1: Consider a single source broadcasting blocks of k packets to n users in a rateless transmis- 
sion. Each packet is a vector of length m over a finite field F^. The channel between the source and each 
of the users is a time invariant erasure channel with probability of erasure p. In each time slot, the source 
broadcasts a random linear combination of k packets. Using random linear coding arguments introduced 
by Ho et al. Q, for a large enough field size d, it is sufficient for the users to receive approximately k 
coded packets to be able to decode the block. 

Let r[t] represent the probability that any given user receives at least k packets in t > k time slots. 
Then, r[t] = $^*=fc (/) (•!■ —pfp^ 1 . Here r[t] represents the fraction of users that have successfully 
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decoded k packets by time t. 

To compare the behavior of r[t] as a function of t for different values of k, we define a normalized time 
variable, s = Accordingly, we define r'[s] = r[ks + k], which can be interpreted as the fraction 
of users that have successfully decoded a single packet in a block of k packets by s time slots. The 
comparison of r'[s] for different k allows us to see, in a normalized time scale, the fraction of users 
that can decode an equivalent of a single packet from a batch of k. We numerically evaluate r'[s] as a 
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Fig. 1. Fraction of users that have successfully decoded a single packet in a block of k packets in s time slots, r'[s] as a 
function of s for p = 1/2 

function of s for different values of k as shown in the Figure Q] for the case where p = 0.5. It can be 
seen from the graph that for k = 10, a large fraction of users are served within a short duration and then 
the source takes a relatively longer time to serve the remaining small fraction of users towards the end 
of the transmission of the current block of k packets. On increasing k to k = 100, the graph becomes 
sharper indicating that the source serves a larger fraction of users in a shorter duration and takes lesser 
time to serve a smaller fraction of users towards the end of the transmission. 

Ideally, we would like all the users to complete decoding together for an increase in throughput. This 
can be achieved by increasing k indefinitely as observed from Figure [TJ However, this causes the decoding 
delay to increase indefinitely as well. Hence, it is important to understand the throughput-delay tradeoff 
as k scales as a function of n. 

III. System Model 

In this work, we study the basic wireless broadcast scenario depicted in Figure [2] that models the 
characteristics of cellular or satellite systems and serves as the fundamental building block for more 
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general networks. In particular, we consider a single source node, S, broadcasting an infinite backlog of 



HOL coding block containing K packets 




Fig. 2. A single source broadcasting to n users over erasure channels with probability of erasure, p in each time slot. 

common information to n users over independent time-varying erasure channels. The data is encapsulated 
into packets, each represented as a vector of length m over a finite field F^. We assume a time-slotted 
operation of the system with Ci[t] G {0, 1} denoting the state of User i's channel in slot t. A single 
packet may be broadcast in each time slot by the source and the transmission to the i th user is successful 
only if Ci[t] = 1. Let C[t] = [Ci[t], . . . ,C n [t]] be the n— dimensional vector of channel states of all 
users in slot t. We refer to C[t] as the channel state of our system. C[t] G S n where S = {0, 1}. For 
simplicity, we assume that all channels are independent and identically distributed. We assume a time 
correlated erasure channel model between the source and each of the users that is defined next. 

Definition 1 (Time Correlated Erasure Channel Model): In each time-slot, user i's channel is in one 
of two states as shown in the Figure [3] When the channel is in the ON state, Cj [t] = 1 with probability 1 
while in OFF state, Ci[t] = with probability 1. The state of the channel evolves as a Markov chain with 
transition probabilities a and (3. Let p denote the steady state probability of an erasure over the channel. 
Then p = ^rg. By setting a + /3 = 1, we can obtain, as a special case, the time invariant erasure channel 
model, where C{ [t] = 1 with probability 1 — p and d [t] = with probability p. 

We consider the class of block coding strategies employed by the source, where data is transferred in 
blocks of k packets. Specifically, the source can start transmitting the next block only if the previous 
block is successfully transferred to all n receivers. Moreover, we focus on the Random Linear Network 
Coding (RLNC) strategy that is defined next. 

Definition 2 (Random Linear Network Coding (RLNC)): In each time slot, the source transmits a ran- 
dom linear combination of the k packets in the Head-of-line (HOL) coding block (see Figure [2]) with 
coefficients of combination chosen uniformly at random from the field F^. In what follows, we refer to 
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Fig. 3. A time correlated erasure channel model. When the channel is in state ON the transmission is successful with probability 
1 while in state OFF the transmission fails with probability 1. 



k as the coding window (or block) size of RLNC. 

Using random linear coding arguments [7 ], the probability that the original k packets can be decoded from 
any k encoded packets formed using the RLNC scheme is equal to ^1 — -jj^ ( 1 ~~ (jfc-i ^ ' ' ' ~ ~d 
which is bounded below by ^1 — ■ Therefore, for a large enough field size d, it is sufficient for the 
users to receive approximately k coded packets to be able to decode the block. Hence, we assume that, 
under the RLNC scheme, the source continues to transmit encoded packets of the current block until 
each user successfully receives k linear combinations. 

It has been shown in |8 ] that random linear network coding is capacity achieving for multicast connections 
in an unreliable network setting. That is, for k sufficiently large, under the coding scheme defined in 
Definition [2j the (broadcast) capacity of our system is (1 — p). Next, we define the two metrics of interest 
in our analysis, namely throughput and decoding-delay. 

Definition 3 ((Broadcast) Throughput): We let R[t] denote the number of packets transmitted by the 
source in a total of t slots. Then, the (broadcast) throughput for a given n and k under RLNC scheme, 
denoted as r](n,k), is the long-term average number of successfully transferred data packets to all n 
users. Hence, we have 

R[t] 

rj(n,k) = lim (1) 

t— >oo t 

Definition 4 (Decoding-Delay): We let denote the number of time slots it takes for the i th user 
to decode the j th block of k packets under the RLNC scheme. The decoding delay of the j block, for 
a given n and k, under the RLNC scheme, denoted as is the time required to transmit all packets 
of the head-of-line (HOL) block to all the users. Hence, we have 

U {j) = max y. (j) (2) 

l<i<n 1 

Recall that the source transmits linear combinations of the current block until each user successfully 
receives k combinations. Then Y^' is the time it takes for k successful transmissions on User i's channel 



s 



and can be written as the sum of durations between each of the k successful transmissions. Hence, 



h=l 



where xjjj is the duration between the (h — l) th and h th successful transmission on User i's channel 
for the j th coding block. 

To understand the importance of scaling k as a function of n to guarantee a non-vanishing throughput, 
consider a simple time invariant channel model, i.e, let Ci[t] be a Bernoulli random variable with p being 
the probability that Ci[t] = in any given time slot t. Owing to the block transmission structure together 
with the independence of the channel states across time, the decoding delay of the block transmissions 
{U^}j>\ are independent and identically distributed. This allows us to model the RLNC operation as 
a renewal process with renewals at the start of each coding block formation and {U^'}j>i being the 
sequence of inter-renewal intervals. Let the random variable U be identically distributed as {U^'}j>i, 
and, E[C/^] = E[f7] for all j. Now, by defining a constant reward of k acquired in each renewal interval, 
we can utilize the main result from renewal theory (9] to write: 

, n ,. R[t] k 

Under the time invariant erasure channel model, it is easy to see that for each i and j, Y> is a negative 
binomial random variable of order k and success probability (1 — p). The exact expression for K[U) is 
as follows, 



i=oo 



E[U]=k+J2 



t=k 



r — 1 



k-1 



p (T-k) qk 



(5) 



\ T=k 

where ( s ) gives the number of size m combinations of s elements and q = (1 — p). 

It can be seen from the above expression that, when A; is a constant independent of n, the mean 
decoding delay E[C7] increases with n. Thus, for any fixed k, rj{n, k) in (0]) goes to zero as n approaches 
oo. However, the exact expression for K[U] is difficult to simplify further. 

In the next section, we provide throughput and delay analysis of RLNC for the time-correlated channel 
model. By identifying a suitable renewal process with reward R(t), we express throughput in terms of 
E[UW] and k using the main result in renewal theory for renewal-reward processes [9]. However, explicit 
characterization of E[[/w] is difficult. Instead, we derive upper and lower bounds on EfCA 7 ')] which enable 
us to understand the scaling of A; as a function of n to guarantee a non-vanishing throughput. Further, 
we compute the mean and variance of when k scales faster than ln(n). 
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IV. Analysis of RLNC for Time Correlated Erasure Channels 

Under the time correlated channel model (when a + /3 7^ 1 in Definition [D, the RLNC operation is not 
a renewal process after each coding block transmission and (0]) is not valid. Hence, in order to express the 
broadcast throughput in terms of the decoding delay, we model the RLNC operation over time correlated 
channels as a semi-Markov process and identify a suitable renewal-reward process. By deriving upper 
and lower bounds on the decoding delay and hence, the throughput, we show that a phase transition in 
the throughput performance of our system occurs at the block-length scaling rate of k = B(ln(n)). This 
result is summarized in the following theorem: 

Theorem 1: Under the RLNC scheme, the broadcast throughput 77(77,, k) of our system, when a, /3 > 0, 
and a + f3 ^ 2, can be characterized as follows, 

a) If k = o(ln(n)), then 

lim 77(71, k) = 0. (6) 

71— >CO 

b) If k = w(ln(n)), then 

lim 77(11, k) = 1 — p. (7) 

n— >oo 

c) Furthermore, if k = 0(m(n)), then 

lim inf 77(71, k) > r(l —p), (8) 

n— >oo 

where < r < 1 is given in (|2TI ). 
In what follows, we develop the mathematical model required to prove the above result. We also 
characterize the asymptotic distribution of the decoding delay for the scaling regime of k = w(ln(n)) 
and derive approximate expressions for the mean and variance of the decoding delay using extreme value 
theory. These expressions reveal the effect of channel correlation on the throughput and delay performance 
of RLNC. 

A. Throughput and Delay Analysis 

Let the random variable G S n denote the channel state of the system at the end of the j th block 
transmission. Let T^> G Z + denote the time slot of the completion of the j th block transmission. Set 
T (o) _ Then? the stoc hastic process {E[t],t > 0}, where, for each j > 0, E[t] = E (j) for t in the 
interval T^> < t < T^ +1 \ is a semi-Markov process where the state transitions occur at the end of each 
block transmission. Let E^ ) = 1, i.e, the process starts with all the user channels being ON. We have 
that = — T^~ x \j > 1 is the time between successive state transitions. Note that U^' is the 
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decoding delay of the j th block transmission and is a non-negative random variable that depends only on 
the current state and the next state of the Markov chain. The embedded Markov chain {E^\ j > 0} is 
an irreducible finite state (2 n states) Markov chain when and f3 ^ 0. Also, the embedded Markov 

chain is aperiodic when a^lor/J^lI^] Let the random variable G S n denote the channel state 
of the system at the beginning of j th block transmission. 

The time epochs of return to the channel state s = 1 in the semi-Markov process form a renewal 
process. The source completes transmission of a random number of blocks in each inter-renewal interval. 

Let the random variable W denote the length of one such inter-renewal interval. Let M denote the random 

M 

number of block transmissions completed in this interval. Therefore, we have the identity W = . 

3=1 

Now, M represents the number of transitions in the embedded Markov chain between consecutive visits 
to state 1. Since the embedded Markov chain is an irreducible aperiodic finite state Markov chain, we 
have that the steady state probabilities of each state s in the embedded chain, 7r s is strictly positive. Hence 

E[M] = i<oo@. 

Now, since E[M] < oo, by defining a reward of 1 for each packet transmitted in a renewal interval, 
we can utilize the main result from renewal theory (9j to write: 

R(t) Mi [Ml 

rj^ k )= lim _Li = IJ (9 ) 

i-s-oo t K[W\ 

The random variable W depends on the channel state of the system at the beginning of each of the M 
block transmissions. Also, owing to the time correlated channel process, the random variables {U^'}jL 1 
are not i.i.d. An explicit characterization of E[W] in terms of \UV'} is, therefore, difficult. Instead, we 
obtain upper and lower bounds on E[W]. Recall that j > 1 is the decoding delay of the j th block. 
Now using (fS]) and ([3]), we can bound U^' as follows, 

max X$ < U {j) < max X$ + max (x^ + ... + X?2\ a.s. (10) 

l<i<n u l<i<n 11 l<i<n \ lZ %k ) 

Note that, for all h > 2, X>0 is the number of time-slots between the (h — l) th and h th successful 
transmission given that the channel state of User i in the time slot of reception of (h — l) th linear 
combination is ON. Therefore, the random variables {X$} i = 1, . . . , n, V? > 1, V/i > 2 are i.i.d. We 
also have that, for all j > 1 and h > 2, X^J is independent of E (1) , . . . ,E ( - J_1) . Also, note that X$ , 
j > 1 is not independent of E^ -1 ^. 



3 The cases when 1) a = = 1 and 2) a — or /3 — are trivial since the channel becomes deterministic within a finite 
expected number of time-slots. 
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At 



Now, the identity W = U^' along with (fTOb enables us to bound E[W] as follows, 

3=1 



E[W] > E 



E[W] < E 



Emax X 
Ki<ra ' 



3=1 
' M 



(3) 
il 



3=1 



max X 

Ki<n 



(3) 



+ E 



3 = 1 



(3) 



EfisOtf + ■■■ + *} 



^(3) 



(11) 



(12) 



To prove TheoremQ] we further bound the lower and upper bounds on K[W] in (fTTT ) and (fT2l ). respectively. 
First, we state a useful lemma iTTOl : 

Lemma 1: (From |[T0l , pg. 7) Let Xii * = 1; 2, . . . , n be identically distributed random variables. Then 
for any r > 0, we have that 



[ max Xi]<- \ Mn) + ln(E[e rXl ])) 

Ki<n T \ 



(13) 



Note that the above bound is useful only when r lies within the radius of convergence of the moment 
generating function (m.g.f) of xi- 

Next, we prove a few lemmas that will be useful in the proof of Theorem [TJ Recall that B^') e S n 
denotes the channel state of the system at beginning of j th block transmission. For any fixed n, \S n \ is 



finite. Hence, maxE 

S65" 



max Xl J h '\B 



Ki<n 



exists and can be bounded as follows, 



Lemma 2: There exist constants fi^ and tq such that, for all j and h, 



maxE 



max Xy h '\B 



0')|r(3) 



Ki<n 



< — (ln(n) + no) . 

TO 



Also, E 



(7] 

max Xf, 

Ki<n in 



< 00 for any j, h G {1, . . . , k}. 



Proof: Fix a r < Miqg)- Then E 



(7) 

exists and is finite. Let R- be the channel state of 



the i th user at the beginning of the j tn block transmission. Let /jq = max E 



;th 



ce{o,i} 

since the channel state of any user is independent of other users' channels, E 
all i, j, and, h G {1, . . . , k}. 



E 



max X y ,'\B 

in 1 



tf)i R (j) 



KKn 



< E 



— In f max e r ° X ^ ) |B (j) = s 

TQ \l<i<™ 



e roX«' + ___ + e r X^| B (3) =s 



- 7 ln v 

< - (ln(n) + ln^o)) 

TQ 



Now, 
< /io for 
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M 



(i). 



In the next lemma, we find an upper bound on the expected value of max X;{ 

— ' l<j<n 
3=1 " " 

Lemma 3: Let \x = — (ln(n) + ln(/uo)) where to and //q are defined in the previous lemma. Then, 



E 



M 



(3) 



x ii 

*■ ' 1<1<TI 



3=1 



< E[Af]/i. 



(14) 



Proof: Let t£$) = max XjjP. Recall that E^ G S n is the channel state of the system at the end of 

l<j<71 

the j th block transmission. Let Tj be the smallest a— field of events containing the a— fields generated 



by the random variables {X£ , . . . , X£ , » = l,...,n,h = l,...,k} and EW,...^^. Let Z = 0. 

m 

Define = - where = E^l^-i]. Also, 

3=1 

ft) = E[^')|jr i _ 1 ] 

= 5^E[C W |B^ = s, ^_i]P(B^ = s^-i) 
= J]E[^)|bO') = s]P(B^ = s|J}_i) 

< Ml 

where the last inequality follows from Lemma [2] Then E[£^")] = Ef/i^)] < [i. We have that {Z^} is a 
martingale by definition. Also, M is a stopping time with respect to the filtration {Fj} with E[M] < oo 
and 

E [| Z (m+l) _ ^M||JT m ] = E [|e (m+1) - ^ (m+1) ||J m ] 

<E[e (m+1) |j m ] + /i (m+1) 

= 2/i( m+1 ) < 2p < oo. 

Hence, using the optional stopping theorem (Theorem 7.2.2 in |9|), we have E[Z( M )] = E[ZW] = 0. 
Therefore, 

r M 



= E[^] = E 



> E 



M 

E« w 

' M 

Ef u > 

3=1 



E 



E« 

_i=i 
E [Mil] . 



U) 



Thus, we have E 



M 

E« M 

3=1 



< E[M]// as required. 
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(i) 

In the next lemma, we find a lower bound on the expected value of max X;i ■ 

l<i<n 

Lemma 4: Let A = — ln(l — j3). Then, for all j, we have 



E 



max X; 



(?) 



KKn 



?2 



a 



> 1 + - [ln(n) + 7] . 
A 



(15) 



Proof: We have that P(X$ = 1) = 1 - a and P(X$ = u) = a/3(l - /3) n " 2 , for all u > 2. Hence, 

u 

PiX^ < u) = (1 - a) + J2 «/?(! " P) v - 2 



l-a(l-/3) 



u=2 
U-l 



Therefore, P( max A^' } < u) = (1 - a(l - /3) M " i ) n Vu > 1 and E 

KKn 



(j) 

max A/, 

KKn iz 



is given as 



E 



max X$ 

Ki<n lA 



l + J^U-h-ae-^-V 



u=l 

coo 



1 - (1 - ae- Xx ) n ) dx 



> 1 + 



By setting w = 1 — ae € [0, 1) in the above integral, we have the following inequalities, 



E 


max X)t 






l<i<n lA 





1 — w r 



dw 



a A(I-IO) 

n-1 

w z dw 

z = 

w z dw 



. rl n-1 

1+ A , 



- 1 + - V 

n-1 n-1 

>i + irJ_.Irl^ 

A^z + l A ^ z + 1 



n-1 



o 



> 1 + - [ln(n) + 7] 
A 



Next, we prove Theorem Q] using the above lemmas and Wald's first equation to bound the lower and 
upper bounds in (TTTb and (fT2l) . respectively. 

Proof: (Proof of the Theorem [D 
Proof of (a) : First, we prove that k = o(ln(n)) leads to a vanishing throughput. As mentioned in 
Lemma [3j M is a stopping time with respect to the filtration {Pj}. We have that { max X^)}j>\ are 

Ki<n — 
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identically distributed and have a finite mean. Also, max Xf> is independent of Tj-\. Thus, using 

l<i<n 

Wald's first equation and Lemma 01 we can lower bound E[W] (cf. CCD) as follows, 



E[W] > E[M]E 



max X 

Ki<n 



(1) 
(2 



n 



> E[M] ^1 + - [ln(n) + 7] J . 
Consequently, using (O, we can upper-bound the throughput as follows, 

m\M}\ 



rj(n, k) < 



E[M] (A + a[ln(n) + 7]) 
kX 



(16) 
(17) 



A + a [ln(n) + 7] 

Clearly, when k = o(ln(ra)), the upper bound, and hence, rj(n, k) decreases to zero, as n — > 00. 
Proof of (b) : Next, we consider the case when k = w(ln(n)). We let k = f(n) ln(n) + 1 for some 
function fin) > such that lim f(n) = 00. Let E (i) = X$ + . . 

n— »oo 

and variance of Y . Then, 



X$ . Let ji and a 2 be the mean 



A = (*-i)(i + |) = (*-i) r ^, 

„ 2 (fc-l)q[2-(a + /3)] 

* = ^ • 

As mentioned in Lemma [3j M is a stopping time with respect to the filtration {J-j}. We have that 
{ max E h>i are identically distributed and have a finite mean. Also, max E ' is independent of 

l<i<n ~ " l<i<n 

Tj-i. Using Wald's first equation and Lemma |3l we can rewrite (fT2l) as follows, 



E[W] < E 



M 

E 

3=1 



max X, 

Ki<n 



a 



+ E[M]/i + <TE 



M 

E 

i=i 



E 



(j) 



max • 

Ki<n 



a 



< E[M] Eu + jx + o-E 



max ■ 

l<j<n 



E (i) - A 



(18) 



Next, we upper bound E 



E 



Cf) 



max • 

Ki<n 



by appealing to Lemma [T] Choose r = 6yQn(n) where b is 



constant such that < r < STn(j3-j). For such a t, 

-r£\ / exp (J) [1 - a (1 - exp (f )) - exp (j) (1 - /?)] \ *~ 1 



E 



exp r- 



E (i) -£ 



(J 



exp 



(7 



l-exp(f) (1-/3) 



r(n) 



/(n) ln(n) 
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where r(n) is defined as follows 



r(n) 



exp 



//(") 



1 — a I 1 — exp 



exp 



(1-/3) 



with di 



-b\fa 



and c?2 



1 — exp 

b/3 



'fin) 



(1-/3) 



y/2-a-fi ^a{2 - a - (3) 



Now, using Lemma [TJ we obtain the following upper bound, 



E 



max 

l<j<n 0" 



< - ln(n) + In I E 



exp t 



= -VlnR(l + /(n)ln(r(n))) 
Using (fT8l) and ( fT9l ), we can bound E[W] and rj(n, k) as follows, 



(19) 



E[W] < E[M) 
r](n, k) > 



<f>(n) 



where (j>(n) is as follows, 
fl-p 



<p(n) 



to 



(ln(n) + ln(/i )) + f(n) ln(n) + 



1 — p 

(l-p)(/(n)ln(n) + l) 
0(n) 

(l-p)v7NV«[2-(« + /3)] 



6/3 



(20) 



ln(n)(l + /(n) In (r(n))) 



Since, f(n) — >■ oo, as n — > oo, we have r(n) — >■ 1 as n — > oo. Let = h ^ + —- Then, 
n->oo /(n) ln(n) + 1 



1 = (I3 lim v/ /(n) In (r(n)) 
= (f 3 lim i//(n)(l - r(n)) In (1 - (1 - r(n))) Tr7 M 



-d 3 lim y/f(n)(l - r(n)) 



(i) 



-da 



exp 



0. 



B n— >oo 



lim V?R 1 - (1 - P) exp 



n) 



1 — a 1 — exp 



(1-/3) exp 



v7( 



(ii) 



1 

where (i) follows by noting that (1 — (1 — r{n))) !-»■(«) — >■ e as n — > 00 and (ii) follows by L' Hospital's 
rule. Hence, from d20l and the fact that rj(n,k) < (1 — p), we see that under the scaling regime of 
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k = w(ln(n)), we have r/(n, k) — > (1 — p) as n — > oo. 

Proof of (c) : Finally, consider the case when k = 0(ln(n)). This is achieved by letting /(n) = b for 
some constant b independent of n in the above analysis. Then, r(n) is a constant independent of n and 



0(n) l-p>Ra+l/ . n 

lim = 1 H — + -e l+61n(r(nj) (21) 

^oo/(n)ln(n) + l To b h ^ h V W V 

> 1 

and from (l20l ). we see that a constant fraction of the capacity is guaranteed, hence proving (c). ■ 

B. Computing the Mean and Variance of Decoding Delay 

In Theorem [1] we showed that it is necessary to scale k at least as ln(n) to guarantee a non-vanishing 
broadcast throughput. Next, we aim to obtain an accurate characterization of the decoding delay, U^' in 
the scaling regime of k = w(ln(n)). In what follows, we drop the superscript (j) and let U denote the 
decoding delay of our system under the RLNC scheme. 

First consider the case of time-invariant channel model. We see that, for each j, Yf' in (0 is the sum 

of i.i.d random variables. By appealing to the Central Limit Theorem, we know that, after suitable 

(i) 

standardization, if' converges in distribution to a standard normal random variable as k — > oo. It is 
also well-known that, the distribution of the maximum of n normal random variables, after suitable 
standardization, converges weakly to the Gumbel distribution ifTTTl as n — > oo. Therefore, we expect the 
decoding delay, U, to converge in distribution to the Gumbel distribution, after suitable standardization, 
as n — > oo. This is shown to be true in lfT2l for the general case of maxima of triangular arrays. We 
summarize the result for our setting in the following theorem, 

Proposition 1: (cf. lfT2l . pg. 961) Let fi(k) = and a 2 (k) = k rrj^p • When k = w(m(n)), we 
have 

lim Pl- 

a{k) 



lim P ( - < a n x + b n \ = exp {-e~ x ) , (22) 



where a n ~ — = and b n ~ a/2 ln(n) as n — > oo. 
V2m(n) 

In general, convergence in distribution does not imply convergence in moments. However, for the case 

of distributions belonging to the domain of attraction of the Gumbel distribution, convergence of moments 

„ ~ U-u(k) 
holds true 02). Let U = Then, Vr > 0, we have, 

a(k) 



lim E 

n— >oo 



U-b n 



r 

■oo 



x r dexp {-e- x ) . (23) 

oo 
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7T 2 



Now, / xdexp (— e x ) = 7 where 7 is the Euler's constant and / x 2 dexp (— e x ) = h 7 2 . 

This enables to write the following approximate expressions for the mean and variance of decoding delay, 
when n is large, 

E[U] ^J^ + ^l ( v / 2 !n7^)+7 ; 1 ), (24) 

1 J i-P i-P v Y v ; V^Ry 

Var[U] W 12(l-rt 2 ln(n) ' (25) 
By using d24l) in (01), we can verify that when /c = u>(ln(n)), the throughput converges to 1 — p as n — > 00. 
In Section [V] using simulations (cf. Figure [6l |7]), we show that the above approximate expressions are 
accurate even for small values of n. 
In lUl, the authors show that, when k is fixed, 

lim P (ll < a n x + & n ) = exp (~e~ x ) , (26) 



where a n = and b n ~ ln(n) as n 



00. 



Now, since a n is a constant, the asymptotic distribution of the decoding delay, U, is concentrated around 
the point b n with a fixed variance for all n and hence, there exists a phase transition with respect to 
the decoding delay, U, such that there exists a threshold on the number of transmissions below which 
the probability that a block of coded packets can be recovered by all the nodes in the network is close 
to zero. On the other hand, if the number of transmissions is slightly greater than the threshold, then 
the probability that every node in the network is able to reconstruct the block quickly approaches one. 
However, when k = w(ln(n)), var[U] in (|25T ) increases with n and the phase transition phenomenon 
observed in [5] becomes less evident. 

For the case of correlated channels, we analyze the decoding delay for the case when the channel state 
of the system is s = 1 in the time slot prior to the beginning of the new block transmission. Let the 
random variable U denote this decoding delay as well. In this case too, the decoding delay for each user 
can be written as a sum of i.i.d random variables. Hence, we have a similar proposition as Proposition Q] 
We can write the approximate expressions for mean decoding delay and variance of decoding delay, when 
n is large, for both the channel models, concisely, as follows, 

* , vr P If 2 ^ (VW)+ 7, (27) 

y(21n(n)) 
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Note that the above expressions are the same as (l24l ) and (l25l ) for the time invariant channel if we 
choose a + j3 = 1 and p = -^s- Although d27l ) is derived under the assumption that the channel state 
of the system before the beginning of the block transmission is s = 1, we show, using simulations (cf. 
Figure [U [9]) in Section |V] it is an accurate upper bound on the decoding delay performance of the actual 
system. 

We can observe an interesting fact from expression d27l) : we see that the mean decoding delay decreases 
as the channel becomes more positively correlated, i.e, 1 < a + /3 — > 2 and increases as the channel 
becomes more negatively correlated, i.e, q + /3<1,q + /3— > 0. The holding time of any state is longer 
when the channel is negatively correlated and hence, the channel remains in a bad state for a longer 
time hurting the decoding delay. On the other hand, the state holding times are shorter for a positively 
correlated channel and state transitions are more frequent leading to a shorter decoding delay. 

V. Numerical Results 

In this section, we provide numerical results to substantiate our analysis of the RLNC scheme, both for 
the case of the time invariant channel model and the time correlated channel model. As a representative 
setup for the case of the time invariant channel model, we let the OFF probability of erasure channels 
p to be 0.1. Note that the broadcast capacity for this choice of p is (1 — p) = 0.9. For the time 
correlated channel model, we choose a = j3 = 0.3. The broadcast capacity for this choice is 0.5. We 
note that the scaling behavior of the throughput and decoding-delay do not change for any other choice 
of channel parameters. Our numerical results are presented under two different scenarios, the first focuses 
on confirming the phase transition of the throughput scaling, and the second focuses on substantiating 
the accuracy of our approximation for the mean decoding delay computed using extreme value theory. 
Study 1) Phase transition: In this study, we explore the phase transition law that is suggested by 
Theorem [TJ To that end, Figure |4] depicts the (broadcast) throughput of RLNC in the actual system 
operation with increasing n for different types of scaling of k for the time invariant erasure channel. We 
see that this result is in perfect agreement with the phase transition law: when k = 150 and therefore scales 
slower than ln(n), we see that the throughput decays towards zero; when k = 501n(n), i.e. k = G(ln(n)), 
the throughput converges to a constant level as suggested by Theorem [Q when k = 101n 2 (n) or n, i.e. 
k = w(ln(n)), the throughput increases toward the broadcast capacity. 

Figure [5] depicts the (broadcast) throughput of RLNC in the actual system operation for the time 
correlated erasure channel model. Once again, we see that this result is in agreement with Theorem Q] 
These two results also reveal that the convergence rate of the performance to the capacity may be 
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Fig. 4. Throughput behavior under different scalings of k 
with n when p = 0.1 



Fig. 5. Throughput behavior under different scalings of k 
with n when a = /3 = 0.3 and capacity is 0.5. 



increased by selecting a faster scaling of k with respect to n. Thus, Study 1 confirms the phase transition 
law suggested by the our analysis. The next study is aimed at studying the accuracy of our approximation 
for the mean decoding delay computed using extreme value theory. 

Study 2 ) Approximate Mean Decoding Delay: Here, we consider two different scalings of k with respect 
to n, and compare the mean decoding-delay of the actual system to the approximate expression obtained 
using extreme value theory. In particular, we study the cases when k = 501n(n) and k = n. 
First, we consider the case of time invariant channel model. Figure [6] depicts mean decoding-delay 
performance when k = 501n(n) of the actual system behavior together with the approximate mean 
obtained using extreme value theory. This demonstrates the accuracy of our approximation even for 
small values of n. We also see that a throughput of approximately 0.85 (see Figure [4]) is achievable with 
this scaling, leading to a decoding delay that scales only logarithmically with the network size. 
In comparison, Figure |7] depicts the mean decoding-delay of the actual system and approximate mean 
when k = n. Again, we observe that our approximation is accurate and applicable to the actual system 
performance, as predicted. In this fast scaling scenario, we also observe that the throughput increases 
towards the capacity of 0.9 (see Figure HJ) instead of converging to a constant level as in the case of 
scaling k as ln(n). Yet, this asymptotic optimality occurs at the cost of linearly increasing decoding-delay 
performance. 

Next, we consider the case of time correlated channel model and compare the actual decoding delay of 
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Fig. 6. Comparing Actual Mean Decoding Delay to 
Approximate Mean for p = 0.1 and k — 501n(n). 



Fig. 7. Comparing Actual Mean Decoding Delay to 
Approximate Mean for p = 0.1 and k — n. 
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Fig. 8. Comparing Actual Mean Decoding Delay to 
Approximate Mean for a = ft = 0.3 and k — 501n(n). 



Fig. 9. Comparing Actual Mean Decoding Delay to 
Approximate Mean for a = /? = 0.3 and k — n. 



the RLNC scheme to the approximate expression in (I27I) . Recall that the approximate decoding delay is 
derived under the assumption that the channel state of the system prior to the beginning of the current 
block transmission is s = 1. Nonetheless, we see, from Figures [8] and |9l that the approximation in (l27l) 
is an accurate characterization of the mean decoding delay of our system. 

Next, we compare the broadcast throughput performance of RLNC for positively and negatively correlated 
channels with the same broadcast capacity of 0.5. From Figure[l0l we see that the throughput performance 
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Fig. 10. Comparing the Rate of Convergence of Throughput when k — n for three cases of Correlated Channels: 1) a = /3 = 0.7, 
2) a = j3 = 0.5, and 3) a = /3 = 0.3. 

is better for the case of positively correlated channels (a = j3 = 0.7) while it is worse for the case of 
negatively correlated channels (a = j3 = 0.3). Overall, these numerical studies collectively confirm the 
accuracy of estimating the mean decoding delay of RLNC using extreme value theory when k scales as 
function of n. 

VI. Extensions 

In this section, we discuss three important extensions of our analysis. First, we discuss how our analysis 
can be extended to LT codes. Next, we comment on the case where the channels are not independent 
across the users. Finally, we comment on the case of asymmetric channels where the channel of user i 
has a steady state erasure probability of p^. 

A. LT Codes 

LT codes Ifl4l are rateless codes designed such that receiving any Vk{$k) = k + O (^/k In 2 (k/ 5k)) 
encoded packets guarantees that the receiver can decode the original k packets with probability (1 — 5k), 
where 5k £ (0, 1). Note that by decreasing 5k, and hence increasing Uk(5k), we can increase the probability 
of successfully decoding at the receiver. Assume that, under the coding strategy of LT codes, the source 
starts transmitting the next block only after all the receivers receive Vk{$k) encoded packets each of the 
current block. The analysis of Section [IV] can be repeated for the case of LT codes by replacing k with 
^fc(^fc) in the derivation of the bounds on the inter-renewal interval K[W] and the computation of mean 
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and variance of the decoding delay. However, in the case of LT codes, the users that receive more than 
Vk(5k) encoded packets have a higher probability of decoding than the other users. Therefore, we have 
to modify our definition of throughput as follows, 

Definition 5 ((Broadcast) Throughput): We let Ri[t] denote the number of packets successfully de- 
coded by the user i in a total of t slots. Then, the (broadcast) throughput for a given n and k obtained using 
the LT coding scheme, denoted as 77(71, fc), is the long-term average number of successfully transferred 
data packets to all n users. Hence, we have 

V (n,k)= lim SXi^M (29) 

The expected reward obtained by the source in one inter-renewal interval of the renewal process defined 

in Section HV-AI is given by — — — where E[i?jl is reward obtained due to user i. Since each user 

n 

can decode the current block with at least a probability of (1 — 8k), we have that, E[Rj] > feE[M](l — 

, fcE[Afl(l-<Sfe) 

Therefore, we have that, 77(77,, k) > — - — = . Also, since the maximum reward that the source can 

E[W] 

obtain from each block transmission is k, we have that 77(71, k) < . If Sp. is such that hx(k/8k) = 

fx \ 

o(k^), then lim = 1. We can now show the following result analogous to Theorem [[} 

fc— s>oo k 

Theorem 2: Using the LT coding scheme, the broadcast throughput 77(71, k) of our system, when a, (3 > 
0, and a + f3 ^ 2, can be characterized as follows, 

a) If k = o(m(n)), then 

lim 77(71, k) = 0. (30) 

n— >oo 

b) If k = u>(m(n)), and if 

i) 8k = 8 for all k, where 8 6 (0, 1), then 

liminf n(n,k) > (1 -p)(l - 8). (31) 

n— >co 

ii) ln(Jfe/* fc ) = o^ 1 / 4 )), and lim 8 k = 0, then 

lim 77(71, k) = (1 — p). (32) 

n— >-oo 

c) Furthermore, if k = 6(ln(n)), and if 

i) 8k = 8 for all k, where 8 G (0, 1), then 

liminf 77(71, k) > r(l — p)(l — 8), (33) 

n— >oo 

ii) \n{k/8 k ) = o(k^ 1 /^), and lim 8 k = 0, then 

k— >oo 

lim inf 77(71, k) > r (1 — p), (34) 
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for some < r < 1. 



Recall that (1 — p) is an upper bound on the throughput achievable by any scheme. A choice of 5k that 



satisfies the conditions 1) \n(k/5k) 



o^C 1 / 4 )) an d 2) lim 5k = is given by 5k 



1 

In(fc) 



Hence, under 



the scaling regime of k = w(ln(n)), this choice of 5k would ensure that the throughput approaches 1—p 

as n — > oo. 



B. Dependent channel model 

Next, consider the case of symmetric time-invariant channels that are dependent across the users. For 
all i, let Ci[t] = with a probability p and C{[t] = 1 with probability Lemma Q] holds true even 

when the random variables x-h * = 1> • • • ; n w& dependent. However, Lemmas [2j [3] and |4] are not valid 
anymore. Nonetheless, we can show that if k = w(ln(n)), then lim r)(n, k) = 1 — p even when the user 
channels are dependent. Next, we prove this sufficient condition to achieve broadcast capacity. 



Once again, we drop the superscript (J) in the following discussion. Let p(k) = j^— and a 2 (k) 



k 



(i-p) : 



-. Define Y^k) 



Note that U can be rewritten as follows, 



U = //(fc) + u{k) max Yi(k) 

Ki<n 



(35) 



Consider the case when k = f(n) ln(n) for some function f(n) > 0. Note that in our analysis we treat 
A; as a continuous function of n. This assumption does not seriously affect the analysis and can easily 



be relaxed. Choose r = b^\a.(n) where b is a constant such that < r < a(k) ln( Such a r lies 
within the radius of convergence of the m.g.f of Yi(k), which is then given by: 



E[e T ^( fc )] = exp 



rA: 



a(k) 
1—p) exp 



exp 



~rp{k) \ 
a(k) J 



1-p 



I-Pexpf^y 



-Wp_ 



1—p exp 



( fe(i-p) 



f(n) ln(n) 



(36) 



Now, using (fL3l with r = 6-^/ln(n), and, (f35l >. and (l36l ). we can bound the expected decoding delay and, 
hence, the throughput as follows: 

g(n) ln(n) 



i—p 

r](n, k) > h(n)(l — p), 



(37) 
(38) 
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where (|38T ) follows from ©, and g(n) and are given as 

( ^ (1 — p) exp 



g(n) = f(n) + 
/(«) 



1 + /(ra)ln 



V 



±vA\ 



1 — p exp 



f b(i-p) 



77(77) 



First, let /(n) = b where b > is a constant independent of n. Then <?(n) and /i(n) are constants 
independent of n with h(n) < 1 (/i(n) can be expressed explicitly in terms of b, b and p). From d38l ), we 
see that the scaling regime of k = 0(ln(n)) guarantees that liminf rj(ri, k) is a non-vanishing fraction 

n— >oo 

of the broadcast capacity (1 — p). Next, let /(n) be such that f(n) — > oo as n — > oo. Then, 



lim 

n->oo /(re) 



lim 



lim ^— v//(n) In 

n— >-oo o 



l + /(n) In 
/ 



V 



1 — p exp 
(1 — p) exp 



( b(i-p) 



6yjA \ 
/(«) / 



lim -^vTfr) 



1 — p exp 
(1 -p)exp 



( t(i-p) 



i 



V 



1 — p exp 



6(1 — p) n^oo 

— lim 

5(1 — p) n->oo 



lim y /(n) 1 — p exp 



f 6(l-p) 

' b(i-p) ' 

.Vp7( 



(1 — p) exp 



77 



v7( 



77 



W"P) /^(l-p)\ -&v^(!-P) 



V7R 



exp 



V vV(«) 



exp 



77 



n 



(iii) 



where (iii) is obtained using V Hospital 1 s rule. 

Hence, we have h(n) — > 1 and 77(77, k) — > (1 — p) as 77 — >■ 00 proving it is sufficient that k = u;(ln(re)) 
to achieve the broadcast capacity (1 — p). 

We can easily see that it is no longer true that if k = o(ln(re)), then lim 7/(77, /c) = 0. For example, 

n— >oo 

consider the trivial case of perfectly correlated user channels. Then the system reduces to the case where 
there is a single user channel and hence, the throughput, 77(77, k), of this system does not depend on 
re. Also, the sufficient condition of k = w(ln(n)) is unnecessary in this trivial case as it is always true 
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that rj(n, k) — > (1 — p) as k — > oo. Although the sufficiency condition appears to be a rather weak 
condition for the trivial example of perfectly correlated channels, it is still useful when the system has 
weak correlations. 

C. Asymmetric channel model 

Consider the case of asymmetric time-invariant channels where Ci[t] = with a probability pi and 
Ci[t] = 1 with probability (1 — pi). The channels of all users are assumed to be independent of each 

other. Suppose that po = suppi exists. Then the broadcast capacity of our system is (1 — po). Let 

i 

E[Ui] and rji (n, k) denote the mean decoding delay and throughput of this system. We can compare this 
system with the symmetric system with erasure probabilities po. Let K[Uo] and rjo(n, k) denote the mean 
decoding delay and throughput of the symmetric system. By comparing the cumulative distributions of 
the decoding delay of the two systems, we can show that E[?7i] < E[C/q]. Consequently, r]o(n,k) < 
r/i(n, fc) < 1 — po- From Theorem [TJ we have that, if k = w(ln(n)), then lim r)o(n,k) = 1 — po and 
hence, lim rji(n,k) = 1 — pq. 

n—>oc 

The case of asymmetric time-correlated channels is technically challenging since it is not easy to obtain 
a closed form expression for the mean decoding delay as in the case of time-invariant channels. 

VII. Conclusion 

We have investigated the throughput and decoding delay performance of RLNC in a wireless broadcast 
setting as the coding window size k scales with the number of receivers n for a time correlated erasure 
channel model. We noted that the broadcast throughput of RLNC vanishes for any fixed k as the system 
size increases. Hence, it is important to understand the scaling of k as a function of n that will guarantee 
a non-vanishing throughput. 

Our analysis revealed a phase transition in the performance of our system, namely, if k increases slower 
than ln(n), the throughput goes to zero as n increases. However, on increasing k faster than m(n), the 
throughput approaches the maximum achievable broadcast throughput of (1 — p). Also, k = 0(ln(n)) 
ensures a constant fraction of the maximum achievable broadcast throughput for the our system. Further, 
we have provided approximate expressions for the mean decoding delay under the scaling regime of k = 
w(ln(n)) using extreme value theory. We have shown through numerical results that our approximation 
is accurate even for small values of n. We have also shown that our analysis can be extended to other 
rateless block coding schemes such as the LT codes. In particular, by choosing 5k such that lim 5k = 0, 
we see that under the scaling regime of k = w(ln(n)), LT codes achieve the broadcast throughput of 
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(1 — p) as n — > oo. Further, we commented on the extension of our analysis for the cases of dependent 
channels across users, and asymmetric channel model. 
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